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Abstract 



In this paper, we study the blowup of the iV-dim Euler or Euler-Poisson equations with 
repulsive forces, in radial symmetry. We provide a novel integration method to show that 
the non-trivial classical solutions (p, V), with compact support in [0, R], where R > is a 
positive constant and in the sense which p(t,r) — and V(t,r) = for r > R, under the 
initial condition 



blow up on or before the finite time T = R 3 /(2Ho) for pressureless fluids or 7 > 1. 

The main contribution of this article provides the blowup results of the Euler (8 — 0) 
or Euler-Poisson (5 = 1) equations with repulsive forces, and with pressure (7 > 1), as the 
previous blowup papers (p] [2], [3] and [4]) cannot handle the systems with the pressure term, 
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1 Introduction 

The isentropic Euler (5 = 0) or Euler-Poisson (S = ±1) equations can be written in the following 
form: 

Pt+V ■ (pu) =0, 

' p[u t + (u • V)u]+VP =pV$, ( 2 ) 
A$(t,x) = Sa(N)p, 

where a(iV) is a constant related to the unit ball in R N : a(l) — 1, a(2) = 27r and a(3) = Air. 
And as usual, p = pit, x) > and u = u(t, x) S R w are the density and the velocity respectively. 
P = P{p) is the pressure function. The 7-law can be applied on the pressure term P(p), i.e. 

P{p)=Kp\ (3) 

which is a common hypothesis. If the parameter is set as K > 0, we call the system with pressure; 
if K = 0, we call it pressureless. The constant 7 = cp/c v > 1, where cp, c v are the specific heats 
per unit mass under constant pressure and constant volume respectively, is the ratio of the specific 
heats, that is, the adiabatic exponent in the equation ([3]). In particular, the fluid is called isother- 
mal if 7 = 1. If K > 0, we call the system with pressure; if K = 0, we call it pressureless. 
In the above systems, the self-gravitational potential field $ = $(t, x) is determined by the density 
p itself, through the Poisson equation ^3. 

When 5 = — 1, the system can model fluids that are self-gravitating , such as gaseous stars. In 
addition, the evolution of the simple cosmology can be modelled by the dust distribution with- 
out pressure term. This describes the stellar systems of collisionless and gravitational n-body 
systems [5] . And the pressureless Euler-Poisson equations can be derived from the Vlasov-Poisson- 
Boltzmann model with the zero mean free path [6]. For TV = 3 and 8 = — 1, the equations ([2]) are 
the classical (non-relativistic) descriptions of a galaxy in astrophysics. See [7] and [5], for details 
about the systems. 

When 5 = 1, the system is the compressible Euler-Poisson equations with repulsive forces. The 
equation ^ 3 is the Poisson equation through which the potential with repulsive forces is deter- 
mined by the density distribution of the electrons. In this case, the system can be viewed as a 
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semiconductor model. See [5] and [TU] for detailed analysis of the system. 
On the other hand, the Poisson equation © 3 can be solved as 



=<5 / G(x -y)p{t,y)dy, 
Jr n 

where G is Green's function for the Poisson equation in the TV-dimensional spaces defined by 



(4) 



G(x) = { 



\x\, N = l; 
log\x\, N = 2; 
- 1 N > 3. 



(5) 



Usually, the Euler-Poisson equations can be rewritten in the scalar form: 



(6) 



For the construction of the analytical solutions for the systems, interested readers should refer 
to pj], [12], [13], [H] and [15]. The results for local existence theories can be found in [16], [17] 
and [18] . The analysis of stabilities for the systems may be referred to [19] , [20] , [21] , p] , [2] , [3] , 
[22], [13], [23], [24], g] and [25]. 

We seek the radial symmetry solutions 



x x 
p(t, x) = pit, r) and u = — Vit, r) —: —V, 

r r 



(7) 



1/2 

with the radius r — ( ^2iLi x i^j 

For the solutions in spherical symmetry, the Poisson equation @3 is transformed to 



^^^rr (t, x) + {N - 1) r N -' z $ r =a (N) 5pr 



N-2, 



„N-l 



4>, 



a(N)6 f , . N , , 
V ' 1 p(t,s)s N - 1 ds. 



(8) 



(9) 



By standard computation, the Euler-Poisson equations in radial symmetry can be written in the 
following form: 



N - 1 

Pt + V Pr + pV r + P V = 0, 



p(V t + VV r ) + P r (p) = p$ r (p) 



(10) 
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Historically, Makino, Ukai and Kawashima initially denned the tame solutions [I] for outside 
the compact of the solutions 

V t + VV r = 0. (11) 

Following this, Makino and Perthame considered the tame solutions for the system with gravita- 
tional forces [2] ■ After that Perthame discovered the blowup results for 3-dimensional pressureless 
system with repulsive forces 3 (8 — 1). In short, all the results above rely on the solutions with 
radial symmetry: 

V t + VV r =^- f p(t, e^ds. (12) 
r Jo 

And the Emden ordinary differential equations were deduced on the boundary point of the solutions 
with compact support: 

D 2 R. 6JVI 

-pp = -^rr, R(0, Ro) = Ro > 0, R{0, R ) = 0, (13) 

where 4? := V and M is the mass of the solutions, along the characteristic curve. They showed 
the blowup results for the C 1 solutions of the system (|10l) . 

Recently, Chae and Tadmor [4 showed the finite time blowup, for the pressureless Euler-Poisson 
equations with attractive forces (<5 = —1), under the initial condition, 

S := {a G R N \ p (a) > 0, n a (a) = 0, V • u(0,a?(0) < 0} ^ <p, (14) 

where fl is the rescaled vorticity matrix (Sl ij) = |(9iUg— 9jUq) with the notation u = (u 1 , u 2 , u N ) 
in their paper and some point Xq. 

They use the analysis of spectral dynamics to show the Racatti differential inequality, 

D div u 1 , , . o , 

< divu) 2 . 15) 

Dt ~ N y ' y ' 

The solution for the inequality (IT5l blows up on or before T = —N/(V ■ u(0, xq(0))- 

However, their method cannot be applied to the system with repulsive forces to obtain the 

similar blowup result. 

On the other hand, in [24] . we have the blowup results if the solutions with compact support 
under the condition, 

2 ( (p\u\ 2 + 2P)dx < M 2 - e, (16) 
Jn{t) 
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where M is the mass of the solution. 

In this article, the alternative approach is adopted to show that there is no global existence of 
C 1 solutions for the system, ([6]) (5 = or 8 = 1), with compact support without the condition 
(fT4| . We notice that the conditions in our result are different from the works of Engerlberg et. al 



Theorem 1 Consider the N -dimensional Euler (5 = 0) or Euler-Poisson equations with repulsive 
forces (5 — 1) 0). The non-trivial classical solutions (p, V), in radial symmetry, with compact 
support in [0, R], where R > is a positive constant (which p(t, r) = and V(t, r) = for r > R) 
and the initial velocity such that: 



Jo 

blow up on or before the finite time T — R 3 /(2Hq), for pressureless fluids (K — 0) or 7 > 1. 

The solutions (p, u) may lose their regularity, for example the velocity function V € C° only or 
the shock waves appear on or before the finite time T. 

2 Integration Method 

In this section, we present the proof of Theorem[T] The technique of the proof was selected simply 
to deduce the partial differential equations to the Racatti equation, to show the blowup result. 
However, we note our integration method is novel to the studies of blowup for this kind of the 
systems. 

Proof. In general, we show that the p(t, x(t; x)) preserves its positive nature as the mass equation 
©1 can be converted to be 




(17) 



-wt +pv - 



(18) 



with the material derivative. 




(19) 



We integrate the equation (fT5|) : 




(20) 
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for po(xo{0, Xq)) > 0, along the characteristic curve. 

We use the momentum equation (|10|) ? with the non-trivial solutions in radial symmetry, po 0, 
to have: 

V t + VV r + K 1P ^ 2 Pr = $ ri 



v t + ^v 2 ) + K 1P -y- 2 p r = ^ r , 

d 1 

rV t + r fr(.2 V2 ) + K 1 r P Pr = r$ '' 
with multiplying r on the both sides. 

We take integration with respect to r, to the above equation, for 7 > 1 or K > 0: 

f R f R d x f R f R 

/ rV t dr+ / r— (-V 2 ) + / K^rp 1 " 2 p r dr = / r$ r dr, 
Jo J dr 2 J ' J 



(21) 
(22) 
(23) 



r dr { 2 V) + 



rVtdr 



for S > 0. 

It follows with integration by part: 



^ dpT- 1 = 



7-1 



r r 



a(N)Sr 



„jV-l 



p(t,s)s N - 1 ds 



f R rV t dr+ [ R rU l -V 2 )+ f ^-dp^ > 0, 
Jo dr 2 J 7-I 



(24) 

dr, (25) 
(26) 



/ rV t dr-- [ V 2 dr+-\RV(t,R) 2 -0-V(t,0) 2 ]- [ J^L p -'- 1 dr+^- \Rp , -\t, R) - • pf^ft, 0)1 > 0. 
Jo 2 7 2 Jo 7-1 7 ~ 1 

(27) 

The above inequality with the boundary compact condition of V(t, i?) = and /d(i, i?) = 0, becomes 



R rV t dr-\ [ R V 2 dr- f ^ 
2 Jo Jo 7-1 



ft- x dr = 0. 



(28) 



As r and t are independent variables and V is C 1 in the domain [0,R] in the assumption of the 
theorem, we may change the differentiation and the integration as the following: 



4- I" r\ dr \ / V- tU- 



7? 



II ^Vdr 2 -- [ R -V 2 dr 2 > 
dt2j 2J 2r 



n 7^1 
R 



n- l dr > 0, 



#7 
7-l' 



9T _1 dr > 0, 



(29) 
(30) 



for 7 > 1 or K = 0. 

For the non-trivial initial condition p > 0j we have the following differential inequality: 

d 1 rR 



dt 2 



1 1 

Vdr 2 1 " 2 ^ 2 



2./o 2r 



V^r J > 0, 



(31) 



Blowup for Euler and EP Equations 



7 



dt 



1 



1 



— / Vdr z > / — V z dr z > — / V z dr z , 



2r 



2R 



dt 



1 



— / Vdr z > — / V z dr 



2R 



By denoting 



H := H(t) = / rVdr = - Vdr 2 , 
Jo 2 J 

and with the Cauchy-Schwarz inequality, 



v ■ ur 



R 

I I V 2 dr 2 



1/2 



f R Vdr< 



R 



R 

I I V 2 dr 2 



n 

o 

1/2 



1/2 



ldr 2 



< ^ / V 2 dr 2 



R 3 ~ 2R 



the inequality (|32j) becomes 



2tf 2 



± H > — I V 2 dr 2 > 
dt ~ 2R J ~ R 3 ' 

d 2H 2 



(32) 
(33) 

(34) 



(35) 

(36) 

(37) 
(38) 

(39) 
(40) 



With the initial condition: Ho = L rV^dr > 0, we can obtain 

H > - R3H \ . (41) 
- 2H t -R 3 y ' 

Therefore, the solutions blow up on or before the finite time T = R 3 / (2H ). 
This completes the proof. ■ 

Remark 2 For controlled experiments in engineering, fluids are kept in a fixed ball solid container 
with a radial R. Therefore, it requires the compact support condition for t > 0. 



p{t,r) = and V(t,r) = 0, 



(42) 



with r > R. This corresponding condition is called no-slip condition (solid boundary condition) 
^ and 

On the other hand, in computing simulations, the systems are usually coupled with the similar 
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boundary conditions for real applications. Therefore, the condition for compact support (no-slip 
condition) is reasonable in modelings. But for free boundary problems, fluids may not be bounded 
by a fixed volume for all time. Therefore, further research is needed to study the corresponding 
result in future works. 

Remark 3 It is still an open question whether or not there exists time-local C 1 -solution with 
compact support for any initial condition with compact support. On the other hand, if the global 
solutions with compact support whose radii expand unboundedly as time tends to infinity, the dis- 
cussion of this paper can offer no information about this case. 

Remark 4 This article has shed new light on situations with the pressure term. In particular, it 
provides the blowup results of the Euler (6 = 0) or Euler-P oisson (b = 1) equations with repulsive 
forces, and with pressure (7 > 1). This is the main contribution of the article, as the previous 
blowup papers (JJlj JEj, and J4%) cannot handle the systems with the pressure term, for C 1 
solutions. A further refinement for the non-radial symmetry is expected in future studies. 
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